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We investigate the electromagnetic transition probabilities for the helium atom embedded in a 
superstrong magnetic field taking into account the finite nuclear mass. We address the regime 
. 7 = 100 — 10 000 a.u. studying several excited states for each symmetry, i.e. for the magnetic 

1) ' quantum numbers 0, — 1, —2, —3, positive and negative z parity and singlet and triplet symmetry. 

The oscillator strengths as a function of the magnetic field, and in particular the influence of the 
finite nuclear mass on the oscillator strengths are shown and analyzed. 
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PACS numbers: 32.60+i, 32.30.-r, 32.70.-n 



I. INTRODUCTION 



Exposing matter to strong and superstrong magnetic fields (which are fields of the order of 10 5 T and above) dra- 
matically changes its properties and yields new and unexpected phenomena. On the microscopic scale, i.e. for atomic 
. t-H . and molecular systems, magnetic forces have a tremendous influence on the electronic structure and quantum dynam- 
ics (USB This is due to the different appearances of the Coulomb and magnetic forces. From a theoretical point 
^ ( — i of view, strong and superstrong magnetic fields are interesting, because the competing forces prevent a perturbative 
treatment of the problem. Therefore it is necessary to develop and apply new nonperturbative techniques. 

Certain astrophysical objects possess strong and superstrong magnetic fields 0,0, 01- Atmospheres of magnetic 
white dwarfs are exposed to fields of the order of 100 - 10 5 T, magnetic fields in the photosphere of neutron stars 



are of the order of 10 5 -10 10 T. For the interpretation of the spectra of these astrophysical objects a wealth of highly 
, accurate atomic and molecular energies, transition wavelengths and transition probabilities are needed. An example 
' for the analysis of astrophysical spectra of magnetized objects using atomic data in strong fields is the white dwarf 
GrW+70°8247, which represents a cornerstone for the understanding of magnetic white dwarfs in general El B, l9l ITo f . 
C^l ■ Highly accurate data are available for hydrogen in strong magnetic fields since more than a decade pj l2tlll| . This 
system is now understood to a very high degree. However beyond hydrogen, there is significant interest in detailed 
data on heavier elements, such as He, Na, Fe and even molecules. Especially helium plays an important role in the 
atmospheres of magnetic white dwarfs and potentially also neutron stars. The electronic structure of the helium atom 



has been considered by several authors during the last decades [T2, [Tj, UJa, [T6|, UJl US UJJ ■ However most of 
. _ ' the corresponding investigations are restricted to a few states or field strengths. Only a few works provide accuracies, 
^ , that are necessary for astrophysical applications. 

Recently detailed investigations of helium in the strong field regime have been performed, providing the community 
with detailed energy levels, transition wavelengths, and transition probabilities, for a den se g rid of field strengths in 
the range of < 7 < 100 a.u. (one atomic units corresponds to 2.3505 x 10 5 T) |2lll22Ll23ll24j . Numerous symmetries 
.£h ' an d many excited states have been addressed. With the resulting large amount of data it was possible to identify 
the absorption edges of the observational spectrum of the magnetic white dwarf GD229 |25|, 1261 l27j , which have been 
$_i ' unexplained for more than 25 years |28l |29| | . 

At this point also the work by Jones et. al [3(j should be mentioned. They applied a released-phase quantum 
Monte Carlo method in order to evaluate bound state energies and dipole-matrix elements for the ground and a few 
excited triplet states. This has been done for a grid of several field strengths 0.08 < 7 < 800 a.u. 

Addressing the superstrong field regime a further challenge is the problem of the finite nuclear mass. The dominant 
energy correction, caused by the finite nuclear mass is for a field of 10 9 T of the same order of magnitude as the 
binding energy itself. This holds even for the energetically lowest states [3lJ • Therefore effects due to the finite nuclear 
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mass have to be taken into account for a correct description of the structure of the atom. Up to date there are no 
detailed studies about the influence of the finite nuclear mass on the transition rates. 

The purpose of the present article is to provide results on the transition probabilities for helium in the superstrong 
field regime. In Sect, ^we review the expressions for the transition matrix elements and analyze the influence of 
the finite nuclear mass. In Sect. IIIII we provide our results and discuss some particular features of the transition 
probabilities as a function of the field strength. Sect. II Vl provides a brief conclusion and an outlook. 



II. ELECTROMAGNETIC TRANSITION PROBABILITIES FOR FINITE NUCLEAR MASS 

A detailed comparison of theoretical and observational spectra requires not only the energies and transition wave- 
lengths, but also the corresponding oscillator strengths. Selection rules of allowed and forbidden transitions are of 
particular importance. Our investigation focuses on the dominant electric dipole transitions. We will shortly review 
the derivation of the corresponding operators since there are modifications due to the presence of the magnetic field 
as well as the finite nuclear mass. 

Our starting point is the pseudo-separated Hamiltonian HE [U 0H1 using relative coordinates {r^} for the 
electrons with respect to the nucleus in atomic units: 



s{K£ +t, ' + 5 B **-') 2 -R} 



2M, 
1 

In - r 2 \ 



(2) 
(3) 



Here is Ma the total mass of the atom, P denotes the pseudo-momentum, and B the magnetic field vector. 
On the other hand, we have the operator H raa < describing the interaction of the system with the electromagnetic 
radiation field A r , neglecting quadratic terms in A r . It is given in relative coordinates by 

Hrad = Y J (^-P + Pi + \ B * r i] A M) + (4) 




A r (r' N ). (5) 

Here r[ denotes the position vector of electron i and r' N the position of the nucleus in the laboratory frame. The 
radiative part of the electromagnetic field A r (r) reads in quantized form (we consider only the creation of photons): 

Ar( r ) = X! N ( k ) a t,\ e k^ exp(ifcr + iu>t) (6) 

a k\ denotes the creation operator for a photon with wave vector k and wavelength A. ek,\ is the polarization vector 
of the photon, whereas N is an amplitude. In the next step, we will integrate over the center of mass coordinate R 
by calculating the matrix element of H ra< i between two eigenfunctions of the pseudo-momentum P (eigenvalues are 
denoted by Ki and Kf), which are given by expressions of the form: 

-Lexp(-iK.R) (7) 

if we assume an integration volume V . 

The dipole approximation, which reads exp(ifc ■ rj) f=a 1, leads us in first order time dependent perturbation theory 
to the following expression for the transition rates: 

^ = 2irJ2[S( E f- E i-")tKi,K f -kX (8) 

|<i|G+|/>| 2 l (9) 
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where 

and < i\, \f >, denote the electronic initial and final states, respectively. In the following we will assume that 
the wavevector k is much smaller than Ki,Kf, which is well-justified in atomic transitions. Thus using e*(-^-p + 

^j^" 2 B x r) =: Q a , we obtain the following expressions for the electronic transitions: 

P% = < f\Q„\i >, 
4°) - ( ? "\ i ^ f\n u 12 



# = ^y^^- (13) 

These expressions represent the dipole-matrix element, the dipole strength and the oscillator strength, respectively in 
the velocity representation. 

On the other hand we have for the expectation value of the commutator 

<i|[H,r]|/> = <i\^ P +^^Bxr\f> (14) 
- (E t - E f ) < i\r\f >, (15) 

where r := r\ + r-2 and p := pi +P2 are symmetrized one-particle operators. Applying the identity of Eqs. (|14fl . (|15|) 

we arrive at the length representation, that reads, 



p% = 2<f\D a \i>, (16) 

,(<7) 

fi 



4? = 4| < f\D a \i > | 2 , (17) 



# = ^^41 (18) 

where e*r =: D a . These above two representations are equivalent. However in case of numerical calculations, the 
two representations yield in general different results. The relative deviation between the two representations is a 
good measure for the convergence of the computational method. Only results that obey certain consistency criteria 
concerning the length to velocity representations of the transition rates are presented. This ensures in particular the 
gauge independence of our results. 

In the following we will assume a vanishing pseudo-momentum K, which is an appropriate approximation in case 
of slow moving atoms. The basic polarization vectors e a are chosen to be parallel and perpendicular to the magnetic 
field vector, indicated as components z and x ± iy. This leads to the following selection rules for the electromagnetic 
transitions of the helium atom in a magnetic field |3(j : 



\M f - Mi\ = 1 and n 2/ n 2 , = 1, (19) 

or 

\M f - Mi\ = and U Zf U Zi = -1 (20) 

and 

S f -Si = and S ZJ - S Zi = 0. (21) 

Here Eq. H19fl . (|2L)fl describe circular and linear polarized transitions, respectively. 

To understand the influence of the finite nuclear mass we rewrite the expression for the oscillator strength in the 
velocity form Eq. (|18fl as 

ftf =2(E f -E l )\<f\D a \ l >l2. (22) 
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The energy factor (Ef — Ei) plays a significant role, as we will see below. One result of Ref.[3l| is, that effects of the 
mass polarization operators are small and therefore in a good approximation results for finite nuclear mass can be 
expressed in terms of results for infinite nuclear mass: 

E{M ,i) « ^£(00,7) - ^ + ^^(00,7) (23) 

Here E(M ,j) denotes the total energy of an eigenstate for the Hamilton operator of the helium atom for nuclear 
mass Mo and a field strength 7. // := (1 — 1/Mq)^ 1 is a reduced mass. First we will concentrate on transitions which 
do not involve tightly bound states. For the corresponding transitions the last, i.e. third term on the right hand side 
of Eq. i|23|) in general cancels in the energy factor of Eq. 1)22(1. As a consequence, there are two generic cases (no tightly 
bound states involved) for the influence of the energy factor (Ef — Ei): In the case of linear polarized transitions, the 
magnetic quantum numbers Mi and Mf are equal and therefore the energy factor is just scaled by the factor 
compared to the results for an infinite nuclear mass. Typically these oscillator strengths are approximately constant 
as a function of the field strength. Note that the factor 1/// deviates for helium about 1CP 4 from 1. In the case 
of circular polarized transitions the second term on the r.h.s. of Eq. l|23[l (— jM/Mq) becomes important, since the 
magnetic quantum numbers Mi and Mf are different. Therefore the linear term 7/-M0 is added to the energy factor, 
which in general causes an increase of the oscillator strengths compared to results for infinite nuclear mass, of the 
form 

/£(M ,7) « ^7/^(00,7) + 5£HP?i(oo,7)| 2 . ( 24 ) 

However, the typical oscillator strengths for circular polarized transitions decrease according to a power law. We 
note that in case of transitions emanating from tightly bound states the third term on the r.h.s. of Eq. \'1'M becomes 
important and in general modifies the pattern for linear polarized transitions. This is essentially due to the fact that the 
energies of magnetically tightly bound states exhibit an inherently different field dependence than the corresponding 
quantity of non tightly bound states. The above discussed behavior will be observed when discussing our results of 
oscillator strengths in Sect lllll 

Some comments on our computational approach are in order. The calculations are performed using an anisotropic 
Gaussian basis set, which was put forward by Schmelcher and Cederbaum |37| . and which has been successfully applied 
to several atoms, ions and molecules [2lt l22l l23t l24l l3ll 133. l39l lifl |41 | . The corresponding basis functions have been 
optimized for each field strength to solve the one particle problems, i.e. H and Hc + in a magnetic field. We refer the 
reader to Ref.|31| for more details. It has been shown, that this approach yields accurate energies and in particular 
oscillator strength for helium, by comparing with the corresponding data in the literature 0, [23, E3- El • 



III. RESULTS 



In this section, we present and discuss our results on the oscillator strengths of electric dipole transitions of helium 
in the superstrong field regime. In order to label the states, we use the standard spectroscopic notation n 2S+1 M n * . 
Here 2S + 1 indicates the spin multiplicity, M is the magnetic quantum number, H z the z parity, and n the degree 
of excitation in the corresponding symmetry subspace. The accuracy for the reported oscillator strengths is between 
1CP 4 and a few times 10 -2 . 

As discussed in Ref.|31| the number of bound states of the helium atom in superstrong magnetic fields becomes 
finite, i.e. the spectrum terminates, if the effects of the finite nuclear mass are taken into account. Therefore only a 
finite, usually small number of transitions "survive" in the superstrong field regime. On the other hand the ionization 
threshold (He — > He + + e~) is up to date not known exactly due to missing detailed investigations on the moving 
He + ion in a magnetic field. The exact field strength for which a certain state becomes unbound is therefore unknown. 
Since our basis functions cannot properly describe the electronic continuum we report here only on transitions that 
are known to be energetically well-separated enough from the continuum. 

The typical features of oscillator strengths of linear polarized transition, discussed in Sect[H] can be clearly seen 
in Fig. n (a). The oscillator strengths for several transitions stay constant, or change much less than one order of 
magnitude. On the other hand transitions emanating from the tightly bound state 1 1 + can be identified by their 
power law behavior. A completely different pattern belongs to the transitions 3 1 + — > 1 1 0~,4 1 + — > 1 1 CP,4 1 + — > 
2 1 0",5 1 0+ -> 1 1 CT,5 1 + -> 1 2 (T, and 5 X 0+ -> 3 1 0~ , depicted in Fig.Q](b). For a field strength below a critical field 
strength j c as 50 they decrease (this can not be seen in Fig. [2(b)), above j c they increase. Numerical values for 
transition wavelengths and oscillator strengths for a few of the lowest linear polarized transitions /it 2S+1 + — * j/ 2S+1 CP, 
(j,, v — 1,2 are presented in tabled 
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FIG. 1: The absolute value of the oscillator strength \fji\ of the linear polarized transitions i/ 1 + — > /i 1 as a function of 
the field strength 7. (a) On the left hand side from bottom to top (i/,^)=(5,4), (5,5), (4,3), (4,4), (3,2), (3,3), (2,2), (2, 1), 
(4,5), (1,1), (3,5), (2,4), (1,2), (2.5), (1,3), (1,4), (1,5). (b) Oscillator strengths of a group of transitions belonging to these 
symmetry subspaces showing a different field dependency. 
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TABLE I: Wavelengths A in A and absolute value of the oscillator strength in atomic units for a few of the lowest linear polarized 
transitions ^ 2S+1 + — > v 2S+1 0~ . The transition 1 1 + — > 1 1 0~ is an example for a linear polarized transition involving a tightly 
bound state. For the transition 2 3 + — > 1 3 0~ the oscillator strength shows a field dependence deviating from the typical 
behavior. 

We present in Fig. [3 the oscillator strengths as a function of the field strength for the circular polarized transitions 
of the form z/0 + -> //(-1) + . The typical power law dependence of the oscillator strengths is observed, as described 
in Sect[H] Furthermore for several transitions we obtain //j(7) ~ C7~ A with a similar exponent A, i.e. parallel 
curves on a double logarithmic scale. On the other hand the reader observes that the number of transitions decreases 
with increasing field strength, being a consequence of the finite nuclear mass effects. Transition wavelengths and 
oscillator strengths for the transitions 1 1 0+ -> 1 1 (-1) + , l x 0+ -> 2 1 (-1) + , 1 3 0+ -> 1 3 (-1) + , and 2 3 0+ -> 1 3 (-1) + 
and oscillator strengths for the transition 1 1 + — > 1 1 (— 1) + with finite mass effects excluded are presented in table ITT1 

In Figs. [31 — El the oscillator strengths of linear and circular polarized transitions are shown as a function of their 
wavelengths for different field strengths addressing the symmetry subspaces v 1 + , i^ 3 + , v 1 ^— 1)~, f 3 (— 1)~, v x {— 2) + , 
and v 3 (— 2) + for v = 1, . . . , 5. With the exception of Figs {7\ and 00 the range of wavelength shown is 10 3 A - 10 5 A. 

Let us first discuss the oscillator strengths emanating from the singlet states with zero magnetic quantum number 
and positive z parity (Fig.[3J|. With increasing field strength the transition wavelengths of some transitions decrease, 
whereas it increases for others. E.g. at w 3000 A a gap between two groups of oscillator strengths emerges and widens 
with increasing field strength. The reader should note, that the values of the oscillator strengths correspondingly 
decrease. Similar statements hold also for the spectrum of the triplet transitions shown in Fig. 

The transition spectrum emanating from the states with magnetic quantum number — 1 and negative z parity shows 
a completely different pattern (see Figs. and . The spectra are only reported up to 7 = 2000, since above this 
field strength there are no transitions between bound states including states of the 2S+1 (— 1) _ symmetry. It can be 
observed that at 7 = 100 there are several very dominant transitions between 10 3 A and 10 5 A (up to 6 atomic units). 
The largest of these disappear with increasing field strength and at 7 = 2000 for 1 (— 1)~, and 7 = 1000 for 3 (— 1)~ 
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TABLE II: Wavelengths A in A and absolute values of the oscillator strength in atomic units for the circular polarized transitons 
l x + -» 1 1 (-1) + , 1 1 + -> 2 1 (-1) + , 1 3 0+ -> 1 3 (-1) + , and 2 3 0+ -> 1 3 (-1) + . Furthermore the table includes fixed nucleus 
results for the oscillator strengths of the transition 1 1 + — > 1 1 ( — 1) + . 
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FIG. 2: The absolute value of the oscillator strength \ffi\ of the circular polarized transitions v + — > 1) + as a function of 
the field strength 7. On the left hand side from top to bottom the following transitions are shown (p, v) — (1, 1),(1, 2), (2, 2), 
(2,1), (1,3), (1,4), (3,1), (1,5), (4,1), (3,2), (2,3), (5,1), (3,4), (4,2), (2,4), (4,3), (2,5), (5,2), (5,3). 



respectively, only one transitions with an oscillator strengths of the order of one remains. 

In Figs. [7| and 00 the spectra for transitions emanating from the 2S+1 (—2) + symmetry subspaces, are presented. 
For 7 = 100, 200 a.u. the oscillator strength increase (with a few exceptions) monotonically as a function of the 
wavelength. For wavelengths of the order of 10 3 A, we find only oscillator strengths much smaller than 1, whereas for 
wavelengths in the interval 10 4 — 10 5 A the corresponding quantities are in the range 3 — 5. At 7 = 10000, only one 
transition of the order of 0.002 atomic units for a wavelength of approximately 10 2 A remains. 

IV. BRIEF CONCLUSIONS 

We have applied a full configuration interaction method to the helium atom in the superstrong field regime between 
100-10 000 atomic units. The effects of the finite nuclear mass have been taken into account. In this work we have 
presented results on the oscillator strengths between bound states. The operators, describing the dominating electric 
dipolc transitions in the magnetic field in first order perturbation theory have been derived from first principles. It 
has been shown how the spectrum changes for different symmetries with increasing field strength. Finite nuclear mass 
effects decrease the number of bound state transitions in the superstrong field regime, since many states enter the 
continuum beyond a certain critical field strength. The influence of the finite nuclear mass on the oscillator strength 
has been analyzed. 

For linear polarized transitions that do not involve tightly bound states the corresponding oscillator strengths are 
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FIG. 3: The oscillator strengths of the linear and circular polarized transitions emanating from the singlet states with 
zero magnetic quantum number and positive z parity, i.e. n 1 + ,n = 1,...,5 with their wavelength given in A for 7 = 
100, 200, 1000, 2000, 10000 a.u. from bottom to top. 
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FIG. 4: The oscillator strength of the linear and circular polarized transitions emanating from the triplet states with 
zero magnetic quantum number and positive z parity, i.e. n 3 + ,n = 1,...,5 with their wavelength given in A for 7 = 
100, 200, 1000, 2000, 10000 a.u. from bottom to top. 



approximately ficld-indcpcndcnt and, compared to the results for infinite nuclear mass scaled by a factor involving the 
reduced mass. For linear polarized transitions involving tightly bound states the oscillator strengths obey a power law 
decay //»(7) ~ C-f~ x . A similar statement holds for the circular polarized transitions. Particular linear and circular 
polarized transitions do not belong to these two cases: they show a different strongly nonlinear dependence on the 
field strength. 

Our results could be of relevance to the interpretation of spectra of neutron stars. For the future an investigation of 
the continuum would be very promising, particular since the discrete spectrum becomes very sparse in a superstrong 
field. The inclusion of motional electric fields into our study would also be very desirable. The latter requires however 
a major theoretical effort. 
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FIG. 5: Same as in Fig. [3] but transitions emanating from the singlet states with n 1 ( — l)~,n = 1,...,5 for 7 = 
100, 200, 500, 1000, 2000 a.u. from bottom to top. 
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FIG. 6: Same as in Fig. but transitions emanating from the triplet states with magnetic quantum number —1 and negative 
z parity for 7 = 100, 200, 500, 1000, 2000 a.u. from bottom to top. 



V. ACKNOWLEDGMENTS 



The Deutsche Forschungsgemeinschaft is gratefully acknowledged for financial support. Discussions with S. Jordan, 
D. Wickramasinghe and J. Liebert are gratefully acknowledged. 



[1] H. Friedrich and D. Wintgen, Phys. Rep. 183, 37 (1989). 

[2] H. Ruder, G. Wunner, H. Herold, and F. Geyer, Atoms in strong magnetic fields (Springer Verlag, Berlin, 1994). 

[3] P. Schmelcher and W. Schweizer, eds., Atoms and Molecules in Strong External Fields (Plenum Press, New York, 1998). 

[4] J. Angel, Ann. Rev. Astron. Astrophys. 16, 487 (1978). 



9 



0.002 
0.001 


0.02 


= 0.2 


2 


4 

2 





7=ioqoo 



7=2000 



y=1000 



y=200 



7=100 



10' 



10 10 
Wavelength [A] 



10 J 



FIG. 7: Same as in Fig. [3] but transitions emanating from the singlet states with magnetic quantum number —2 and positive 
z parity, i.e. n 1 (-2) + , n = 1, . . . , 5 with their wavelength given in A for 7 = 100, 200, 1000, 2000, 10000 a.u. from bottom to 
top. Note the different length scales on the oscillator strength axis, for different field strengths. 



10 ' 



10 



=1 1 

CO 



7=10000 




7=2000 


7=1000 

1 




7=200 

1 






7=100 

1 

'2 '3 


. L4 







Wavelength [A] 



FIG. 8: Same as in Fig. but for transitions of corresponding triplet states, i.e. n 3 (— 2) + ,n 



, 5 for 7 



100, 200, 1000, 2000, 10000 atomic units from bottom to top. Note the different length scales on the oscillator strength axis, for 
different field strengths. 



[5] G. G. Pavlov, Y. A. Shibanov, V. E. Zavlin, and R. D. Meyer, in Proc. NATO ASI C 450, edited by M. A. Alpar, 
U. Kiziloglu, and J. van Paradijs (Kluwer, Dordrecht, 1995), pp. 71-90. 

[6] D. T. Wickramasinghe and L. Ferrario, Pub. Astron. Soc. Pac. 112, 873 (2000). 

[7] J. Angel, J. Liebert, and H. S. Stockmann, Astrophys. J. 292, 260 (1985). 

[8] J. L. Greenstein, R. Henry, and R. F. O'Connel, Astrophys. J. 289, L25 (1985). 

[9] G. Wunner, W. Rosner, H. Herold, and H. Ruder, Astron. Astrophys. 149, 102 (1985). 
[10] D. T. Wickramasinghe and L. Ferrario, Astrophys. J. 327, 222 (1988). 
[11] Y. P. Kravchenko, M. A. Liberman, and B. Johansson, Phys. Rev. A 54, 287 (1996). 
[12] R. O. Mueller, A. Rau, and L. Spruch, Phys. Rev. A 11, 789 (1975). 
[13] J. Virtamo, J. Phys. B 9, 751 (1976). 

[14] P. Proschl, W. Rosner, G. Wunner, and H. Herold, J. Phys. B 15, 1959 (1982). 



10 



M. Vincke and D. Baye, J. Phys. B 22, 2089 (1989). 
M. Braun, W. Schweizer, and H. Herold, Phys. Rev. A 48, 1916 (1993). 

G. Thurner, H. Korbel, M. Braun, H. Herold, H. Ruder, and G. Wunner, J. Phys. B 26, 4719 (1993). 
M. V. Ivanov, J. Phys. B 27, 4513 (1994). 

M. Braun, W. Schweizer, and H. Elster, Phys. Rev. A 57, 3739 (1998). 
J. S. Heyl and L. Hernquist, Phys. Rev. A 58, 3567 (1998). 
W. Becken, P. Schmelcher, and F. Diakonos, J. Phys. B 32, 1557 (1999). 
W. Becken and P. Schmelcher, J. Phys. B 33, 545 (2000). 
W. Becken and P. Schmelcher, Phys. Rev. A 63, 053412 (2001). 
W. Becken and P. Schmelcher, Phys. Rev. A 65, 033416 (2002). 
R. F. Green and J. Liebert, Publ. Astr. Soc. Pac. 93, 105 (1980). 
G. D. Schmidt, W. B. Latter, and C. B. Foltz, Astrophys. J. 350, 768 (1990). 
G. D. Schmidt, R. G. Allen, P. S. Smith, and J. Liebert, Astrophys. J. 463, 320 (1996). 
S. Jordan, P. Schmelcher, W. Becken, and W. Schweizer, Astron. Astrophys. Lett. 336, L33 (1998). 
S. Jordan, P. Schmelcher, and W. Becken, Astron. Astrophys. 376, 614 (2001). 
M. D. Jones, G. Ortiz, and D. M. Ceperley, Phys. Rev. A 59, 2875 (1999). 
O.-A. Al-Hujaj and P. Schmelcher, Phys. Rev. A 67, 023403 (2003). 
J. W. E. Lamb, Phys. Rev. 85, 259 (1952). 

J. E. Avron, I. W. Herbst, and B. Simon, Ann. Phys. 114, 431 (1978). 
B. R. Johnson, J. O. Hirschfelder, and K. H. Yang, Rev. Mod. Phys. 55, 109 (1983). 

P. Schmelcher, L. S. Cederbaum, and U. Kappes, in Conceptual Trends in Quantum Chemistry, edited by E. S. Kryachko 
and J. L. Calais (Kluwer, Dordrecht, 1994), pp. 1—51. 
W. Becken, Ph.D. thesis, Universitat Heidelberg (2000). 
P. Schmelcher and L. S. Cederbaum, Phys. Rev. A 37, 672 (1988). 

U. Kappes and P. Schmelcher, Phys. Rev. A 54, 1313 (1996), 53, 3869 (1996); 51, 4542 (1995). 
T. Detmer, P. Schmelcher, F. K. Diakonos, and L. S. Cederbaum, Phys. Rev. A 56, 1825 (1997). 

T. Detmer, P. Schmelcher, and L. S. Cederbaum, Phys. Rev. A 57, 1767 (1998); 61, 043411 (2000); 64, 023410 (2001); J. 
Chem. Phys., 109 (1998); J. Phys. B 28 (1995). 
[41] O.-A. Al-Hujaj and P. Schmelcher, Phys. Rev. A 61, 063413 (2000). 



